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A new modeling approach for large-eddy simulation (LES) is obtained by combining a 'regu- 
larization principle' with an explicit filter and its inversion. This regularization approach allows a 
systematic derivation of the implied subgrid-model, which resolves the closure problem. The central 
role of the filter in LES is fully restored, i.e., both the interpretation of LES predictions in terms 
of direct simulation results as well as the corresponding subgrid closure are specified by the filter. 
The regularization approach is illustrated with 'Leray-smoothing' of the nonlinear convective terms. 
In turbulent mixing the new, implied subgrid model performs favorably compared to the dynamic 
eddy-viscosity procedure. The model is robust at arbitrarily high Reynolds numbers and correctly 
predicts self-similar turbulent flow development. 
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Accurate modeling and simulation of turbulent flow is a topic of intense ongoing research jjj . Modern strategies 
for turbulent flow are aimed at reducing the dynamical complexity of the underlying system of partial differential 
equations while reliably predicting the primary flow phenomena. In large-eddy simulation (LES) these conflicting 
requirements are expressed by coarsening the description on the one hand and subgrid modeling on the other hand. 
The coarsening is achieved by spatial filtering which externally specifies the physical detail that will ideally be 
retained in the LES solution. Maintaining the dynamical properties of the resolved large scales is approached by 
introducing subgrid modeling to deal with the closure problem that arises from filtering the nonlinear terms. 

In the filtering approach to incompressible flow the specification of the basic convolution filter L is all that is required 
to uniquely define the relation between the unfiltered and filtered flow field as well as the closure problem for the 
so-called turbulent stress-tensor Ty . This situation is in sharp contrast with actual present-day large-eddy modeling 
in which the specification of the subgrid model for Tij as well as the comparison with reference direct numerical 
simulation (DNS) results is performed largely independent of the specific choice of the filter L. 

In this paper we will formulate an alternative approach to large-eddy simulations which completely restores the 
two central roles of the basic filter L, i.e., providing an interpretation of LES predictions in terms of filtered DNS 
results as well as fully specifying all details of the subgrid model. The key elements in this new formulation are a 
'regularization principle', a filter L and its (formal) inverse operator denoted by L^ 1 ||. 

A regularization principle expresses the smoothing of the dynamics of the Navier-Stokes equations through a specific 
proposal for direct alteration of the nonlinear convective terms. This modeling differs significantly from traditional, 
less direct approaches, e.g., involving the introduction of additional eddy- viscosity contributions 0|. The latter are 
clearly of a different physical nature and do not fully do justice to the intricate nonlinear transport structure of the 
filtered Navier-Stokes equations. The regularization principle gives rise to a basic mixed formulation involving both 
the filtered and unfiltered solution. Application of L and L~ 1 then allows to derive an equivalent representation solely 
in terms of the filtered solution. This provides a unique identification of the implied subgrid model without any further 
external (ad hoc) input or mathematical-physical considerations of the closure problem. The regularization modeling 
approach is not only theoretically transparent and elegant, but it also gives rise to accurate LES predictions. In 
particular, we consider the implied subgrid model that arises from Leray's regularization principle A comparison 
between the Leray model and dynamic subgrid modeling (e.g., Q) will be made for turbulent mixing flow, both at 
moderate and at high Reynolds numbers. 

In the filtering approach one assumes any normalized convolution filter L : ui — > Mi where Ui (ui) denotes the 
filtered (unfiltered) component of the velocity field in the Xi direction. Filtering the Navier-Stokes equations yields 

d t ui + dj(TijUi) + dip - -^-djjUk = djUj (1) 
He 

where the turbulent stress tensor Tij = TTiUj — uiUj represents the closure problem and Re denotes the Reynolds 
number. Both the relation between Ui and Ui as well as the properties of Tij are fully specified by L. In actual subgrid 
modeling for LES, the next step is to introduce a subgrid model ra,j(u) to approximate Tij. A variety of subgrid 
models has been proposed to capture dissipative, dispersive or similarity properties of . 

Many subgrid models are arrived at through a physical or mathematical reasoning which is only loosely connected 
to a specific filter L. As an example, the well-known Smagorinsky model [Q is given by mfj — — (CsA) 2 \Sij(u)\Sij(vL) 
where the rate of strain tensor Sij = diUj+djUi and \Sij \ 2 = SijSij /2. The only explicit reference to the filter, made in 
this model, is through the filter-width A. In actual simulations A is specified in terms of the grid-spacing h rather than 
in terms of L. Furthermore, the Smagorinsky constant Cs is determined independent of L, which further reduces any 
principal role for the filter. The situation is comparable for the 'tensor-diffusivity' model = CxD^idkUidkUj , 

with Afc the filter-width in the x^-direction Q. The coefficient Ctd is usually related to the normalized second 
moment (L(x 2 ) — x 2 ) / A 2 of the filter L. For various popular filters such as the top-hat or the Gaussian filter one finds 
Ctd — 1/12, i.e., independent of the actual filter used. The role of the filter is in principle fully explicit in Bardina's 
similarity model mf- = u^Uj — uiUj ||. In actual simulations, however, one frequently adopts a wider explicit filter 
or a filter of a different type, to enhance smoothing properties of this model I]]. Moreover, the model is sometimes 
multiplied by a constant Cb which is specified independently of any presumed filter Jt|. Finally, the successful dynamic 
subgrid modeling requires only the explicit specification of the so-called test-filter |L0| . To retain the central Germano 
identity the test-filter can in principle be chosen independent of L, mainly requiring the specification of the filterwidth 
of the test-filter relative to A. Additional averaging over homogeneous directions, 'clipping' steps to stabilize actual 
simulations, and the fact that the assumed base-models are themselves only loosely connected to L, also make the 
dynamic procedure rather insensitive to the specific assumed filter. 

In contrast to these popular LES models, the regularization approach involves the introduction of a pair (L, L^ 1 ) 
to fully specify the implied subgrid model as well as the interpretation of LES predictions in terms of reference 
DNS results. The selection of any other pair (£, C^ 1 ) directly leads to its corresponding DNS interpretation and 
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the associated subgrid model consistent with the regularization principle. This modeling strategy has a number of 
important benefits, addressing directly the nonlinear convective contributions and requiring no additional 'external' 
information such as model coefficients or the width of the test-filter. The regularization principle allows a transparent 
modeling in which the modeled system of equations can be made to share a number of fundamental properties with the 
Navier-Stokes equations, such as transformation symmetries, Kelvin's circulation theorem, etc.. The implied subgrid 
model is quite simple to implement, with some technical complications arising from the construction of an accurate 
inverse operator L^ 1 . 

To illustrate the approach we consider the intuitively appealing and particularly simple Leray regularization in 
which the convective fluxes are replaced by HjdjUi, i.e., the solution u is convected with a smoothed velocity u. 
Consequently, the nonlinear effects are reduced by an amount governed by the smoothing properties of L. The 
governing equations in the Leray formulation can be written as pi 



djUj — ; dtUi + ujdjUi + dip — Tr^jj u i = (2) 

Re 



Uniqueness and regularity of the solution to these equations have been established rigorously H . The Leray formu- 
lation contains the unfiltered Navier-Stokes equations in the limiting case L — ► Id, e.g., as A — ► (Id denotes the 
identity). The unfiltered solution can readily be eliminated from (|^) by using the inversion operator Uj — L — 
After some calculation (||) can be written in the same way as the LES 'template' (|l|) in which T\j on the right hand 
side is replaced by the asymmetric, filtered similarity-type Leray model mf- given by: 



m% = L 



(v,jL 1 (ui)j — UjUi = UjUi — UjUi (3) 



This model requires the explicit application of both L and L _1 . 

In the sequel we consider invertible numerical quadrature approximating the top-hat filter. In one dimension the 
numerical convolution filtering u = G * u corresponds to kernels 

G(z) = J2^6(z- Zj ) ; N<A/2 (4) 

In particular, we consider three-point filters with ciq = 1 — a, a\ = a_i = a/2 and zq = 0, Z\ = —Z—\ = A/2. Here we 
use a = 1/3 which corresponds to Simpson quadrature of the top- hat filter. In actual simulations the resolved fields 
are known only on a set of grid points {x m }m=o- The application of L~~ x to a general discrete solution {u(x m )} can 
be specified using discrete Fourier transformation as pT| 



where the subgrid resolution r = A/h is assumed to be even. An accurate and efficient inversion can be obtained with 
only a few terms, recovering the original signal to within machine accuracy with n rs 10. At fixed A, variation of the 
subgrid resolution r allows an independent control over flow-smoothing and numerical representation |]l2| . Filtering 
and inversion in three dimensions arises from composing three one dimensional filters. 

To assess the Leray model the turbulent mixing layer is simulated in a volume (? at various Re adopting a fourth 
order accurate spatial discretization and explicit Runge-Kutta time-stepping. We compare predictions with those 
obtained using the dynamic subgrid model, which was shown to be among the most accurate models in a comparative 
study of the same turbulent mixing layer reported in || . 

A first introductory test of the Leray model is obtained by studying instantaneous solutions. As a typical illustration 
of the mixing layer the DNS prediction of the normal velocity U2 is shown in the turbulent regime in Fig. 0(a). We 
used Re = 50 based on the initial momentum thickness and free-stream flow properties. The filtered ui can be 
seen in Fig. |l|(b) establishing a significant smoothing due to the 'Simpson' filter at A = ^/16. The Leray prediction 
(Fig. |l|(c)) appears to capture the main 'character' as well as some of the details of the filtered DNS solution. A 
slight underprediction of the influence of the small scales is, however, apparent. Further visualization showed that the 
instantaneous Leray predictions display much better overall agreement with filtered DNS than the dynamic model, 
which relative to the Leray model significantly overpredicts the smoothing Of course, assessing the quality of LES 
predictions in this way is difficult to quantify and we consider more specific measures next. 

The evolution of a crucial mean-flow property such as the momentum thickness is shown in Fig. |^. The Leray 
results compare significantly better with filtered DNS results than those obtained with the dynamic model on 32 3 
grid-cells. We observe that some of the discrepancies between Leray and filtered DNS results are due to numerical 
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contamination. By increasing the resolution at fixed A, a good impression of t he g rid-independent solution to the 
modeled equations can be inferred using 64 3 - 96 3 grid-cells, i.e., A/h = 4 to 6 [12]. Numerical contamination also 
plays a role in the dynamic model. The grid-independent solution corresponding to the dynamic model appears less 
accurate than the corresponding Leray result. 

A more detailed assessment is obtained from the streamwise kinetic energy spectrum shown in Fig. |^. The dy- 
namic model yields a significant underprediction of the intermediate and smaller retained scales, particularly for the 
approximately grid-independent solution. The Leray predictions are much better. On coarse grids, an overprediction 
of the smaller scales is apparent due to interaction with the spatial discretization method. At proper numerical 
subgrid resolution the situation is considerably improved and the Leray model is seen to capture all scales with high 
accuracy. A slight, systematic underprediction of the smaller scales remains, consistent with the impression obtained 
from Figs. 0(b)- (c). 

A particularly appealing property of Leray modeling is the robustness at very high Reynolds numbers, cf. Fig. |J. 
This is quite unique for a subgrid model without an explicit eddy-viscosity contribution. Although comparison 
with filtered DNS data is impossible here, we observe that the smoothed Leray dynamics is essentially captured as 
r = A/h > 4 j^J. The tail of the spectrum increases with Re, indicating a greater importance of small scale flow 
features. Improved subgrid resolution shows a reduction of these smallest scales, consistent with the reduced numerical 
error. At high Re the spectrum corresponding to the Leray model tends to contain a region with approximately k~ 5 / 3 
behavior, which is absent at Re — 50. Further analysis showed that the solution develops self-similarly at high Re. 

The Leray model was presented to illustrate the new regularization approach for LES. It predicts mean flow 
properties such as momentum thickness very accurately. The model exhibits both positive and negative production of 
turbulent kinetic energy. The computational overhead associated with the Leray model can be much lower than that 
of dynamic (mixed) models, especially if quantities are desired which are rather insensitive to the inversion quality. 
The regularized Leray dynamics shows an appealing robustness at high Re. Further extensions of the regularization 
approach are presently being considered. Of particular interest is the Lagrangian averaged NS— a model [ fi"3[ which 
arises in the Euler-Poincare framework for smoothed flow dynamics. 
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FIG. 1. Normal velocity component 112 at time t = 80, (a): DNS, (b): filtered DNS, (c): Leray on 64 3 ; using a filterwidth 
A = 1/16. The light (dark) isosurfaces correspond to 112 = 0.3 (—0.3). 



FIG. 2. Momentum thickness 9: filtered DNS (o), Leray- model (32 3 : dash-dotted, 64 3 : solid, 96 3 : A), dynamic model (32 3 : 
dashed, 64 3 : dashed with o). A fixed filterwidth of £/l6 was used. 



FIG. 3. Streamwise kinetic energy spectrum E at t — 75: filtered DNS (o), Leray-model (32 3 : dash-dotted, 64 3 : solid, 96 3 : 
A), dynamic model (32 3 : dashed, 64 3 : dashed with o). A fixed filterwidth of £/16 was used. 

FIG. 4. Streamwise kinetic energy spectrum E at t = 75 predicted by the Leray model: Re = 50 (64 3 : dash-dotted, 96 3 : 
dash-dotted, A), Re = 500 (64 3 : dashed, 96 3 : dashed, A), Re = 5000 (64 3 : solid, 96 3 : solid, A). A fixed filterwidth of £/I6 
was used. The dotted line represents fc~ 5 / 3 . 
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FIG. 1. Normal velocity component ui at time t = 80, (a): DNS, (b): 
A = £/16. The light (dark) isosurfaces correspond to U2 = 0.3 (—0.3). 



filtered DNS, (c): Leray on 64 ; using a filterwidth 
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FIG. 3. Streamwise kinetic energy spectrum E at t = 75: filtered DNS (o), Leray-model (32 3 : dash-dotted, 64 3 : solid, 96 3 : 
A), dynamic model (32 3 : dashed, 64 3 : dashed with o). A fixed filterwidth of £/16 was used. 
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FIG. 4. Streamwise kinetic energy spectrum E at t = 75 predicted by the Leray model: Re = 50 (64 3 : dash-dotted, 96 3 : 
dash-dotted, A), Re = 500 (64 3 : dashed, 96 3 : dashed, A), Re = 5000 (64 3 : solid, 96 3 : solid, A). A fixed filterwidth of £/16 
was used. The dotted line represents fc~ 5//3 . 
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